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Abstract 

We analyze forward-backward charge asymmetry of the lepton production 

in rare decays B — » X s l + l~ and B — > K*l + l~, including vector-resonance 

effects. Certain regions of phase space, in which the asymmetry is sensitive 

to individual short-distance coefficients, are pointed out. In particular, we 

suggest a method to test the coupling of the leptonic axial vector current to 

the left-handed quark current experimentally. 
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I. INTRODUCTION 



Rare B decays have been the focus of many experimental and theoretical considerations 
|l||J]. This is due to the amount of information on the standard model that can be extracted 
from these processes. The rare decays proceed through flavor-changing neutral current 
diagrams that are absent at the tree level and thus provide a good probe of the standard 
model at the loop level. For example, recent measurements of photon penguin induced 
processes B — > X s ~f and B —>■ K*^ by the CLEO group [T|] are consistent with the standard 
model predictions. Likewise, the dilepton decay B — > XJ + l~ and corresponding exclusive 
modes will provide tests of the Z boson penguin and W boson box diagrams which take 
over the photonic penguin diagram for large top quark mass. On the other hand, unlike the 
radiative B decays, the dileptonic process are strongly dependent on the top quark mass 
m t through the dominant internal top quark line, so that a comparison between theoretical 
calculations as a function of m t and experimental data may lead to constraints on the 
top quark mass. Such an exercise has been done with the current upper limit of Br (6 — > 
sfi + fi~) < 5.0 x 10~ 5 (90% CL), which gives the bound of ~ 390 GeV on the top quark mass 
0. These rare B decays are also sensitive to quark mixing angles V t d, V ts and Vtb, hence their 
measurements yield valuable information on Cabibbo-Kobayashi-Maskawa(CKM) matrix 
elements and consequently shed some light on CP violation in the standard model. In this 
work we concentrate on the rare dileptonic decay b — > sl + l~ (/ is either an electron or a 
muon.). 

While extensive investigations have been carried out for the invariant mass spectrum 
of dileptons, the authors of Ref. |4[] pointed out that to fully extract short-distance coef- 
ficients both in magnitude and in sign from experiments, one has to consider the angular 
distribution such as forward-backward charge asymmetry of the l + production in the decay 
b — > sl + l~ , either. In addition, this asymmetry can be used to test the chirality of the b — > s 
transition. The process B — > K*l + l~ has been considered for the small recoil case where 
heavy quark symmetries apply [§,0. Recently, the Dalize distribution of the exclusive rare 
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dilepton decays has been derived in terms of various form factors J7| and the sensitivity of it 
to possible new physics has been investigated. In particular, it is shown that measurement 
of the transverse polarization of K* meson yields information on the coupling of the leptonic 
axial vector current to the left-handed quark current. It is the aim of this work to provide a 
comprehensive analysis of asymmetry for both inclusive process B — > X s l + l~ and exclusive 
channel B — > K*l + l~ in the context of the standard model. We attempt to probe kine- 
matic regions in which asymmetries are sensitive to individual short-distance coefficients. 
Thus measurements of these asymmetries provide tests of short-distance physics which is 
sensitive to extensions of the standard model. In this context it is suitable to consider the 
integral asymmetry as we have flexibility to make various phase space cuts and read off the 
corresponding 'partial integrated' asymmetry. 

The forward-backward asymmetry results from the interference of leptons produced by 
the vector current with that by the axial vector current (see Eq. (10)). This is because 
it is a parity violating effect. In the b — > sl + l~, the axial vector current is mainly due 
to the Z boson penguin and W boson box diagrams in which top quarks are involved. 
The vector current arises not only from these two diagrams along with the photon penguin 
one, but also from loops containing instead charm quarks. Extensions of the standard 
model lead to extra diagrams. In models with two Higgs doubles, for example, there is a 
charged scalar particle. The charged scalar coupling to quarks can substitute W boson in 
the Z boson penguin diagram. In addition to these, resonances such as the J/ip and if}', 
created by the neutral four-quark operator bscc, contribute to the vector current via the 
V — > . As shown later on, different sources of the vector current manifest themselves 
in different regions of phase space, whereas the axial vector current is constant overall. 
For the very low q 2 the photonic penguin dominates, while the Z pengiun and W box 
becomes important towards high q 2 . Furthermore, we shall demonstrate in this work that 
the short- distance component of the vector current suffers from cancellation dramatically 
in the middle q 2 region. Actually, the pure top quarks effect vanishes at an invariant mass 
close to the J/ip mass and the analogy can be made with nonresonant contributions near 
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the ip 1 '. Thus the asymmetry over there, roughly speaking, contains to two factors; one is 
from top quarks responsible for the axial vector current, and another is from charm quarks 
for the vector current. The former is short-distance and sensitive to new physics, while the 
latter, we believe, remains unchanged in extensions of the standard model. As it is not 
renormalized under QCD, the effective axial vector vertex can be calculated reliably. On 
the other hand, calculations show that the Wilson coefficient of the bscc operator, denoted 
by 02, is sensitive to the QCD scale parameter and the renormalization point and scheme 
beyond the leading logarithmic approximation ||. But it is hopeful phenomenologically to 
extract this coefficient quite precisely from nonleptonic B decays. With the top mass from 
the Fermilab and the value of determined by the data of the B — > X s J/ip, the standard 
model prediction for forward-backward asymmetry of the middle region is insensitive to 
the dependence on the QCD renormalization parameters. Meanwhile, measurements of the 
asymmetry offer us an opportunity to detect the axial current in the b — > sl + l~ and to 
set constraints on extensions of the standard model. We suggest therefore that one should 
home in on the resonant region for revealing the effect of the axial vector current and short- 
distance physics in general. This is contrast to the decay rate spectrum where one has to 
'cut out' J /if) and if)' resonances to probe the short- distance effect. 

This paper is organized as follows. We will briefly review the effective Hamiltonian for 
the flavor- changing neutral-current-induced transition b — > sl + l~ in the following section. 
Then we proceed into the analysis of forward-backward asymmetry of the inclusive process 
B — > XJ + l~ in Sec. III. The exclusive channel B — » K*l + l~ is also examined in the sec- 
tion, by employing heavy quark symmetry to determine form factors from measurements of 
semileptonic D decays. Sec. IV is devoted to relevant discussions along with our conclusion. 
An Appendix contains some useful integral formulae occurred in the main text. 
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II. EFFECTIVE HAMILTON! AN 



Let us begin with an effective Hamiltonian relevant to flavor- changing one-loop processes 

b^si+r [1, 
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(1) 



with effective vertices, 



r£< B > = A(S) Tm (1 - 75) - im b s 2 w F 2 a, u q»(l + 75)/^ 



The coefficient functions are given by 



A(B,F 2 )= Y, U q A q (B q ,F%), 

q=u,c,t 



(2) 



in which U q = V*V q b is the product of the relevant elements of the CKM matrix. As long as 
the u quark is ignored, one has U c + Ut = from unitarity of the CKM matrix. The heavy 
particles with masses much larger than the physical scale \x ~ m& are integrated out and 
their masses are absorbed into coefficient functions evaluated in two steps pD|HTB|]; first one 
computes these functions at a renormalization point about the heavy mass and then scales 
them down to the order of the b quark mass using the renormalization group equation. In 
the standard model the coefficients A t and B t at the scale equal to the W boson mass, for 



example, take the forms (with x = mf/M^) 
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When the scale goes down, we notice first of all that as the four-fermion operator with the 
axial vector current of leptons has a vanishing anomalous dimension and does not mix with 



other operators because of chirality, the combination A — B retains its value at the higher 
mass scale and thus is independent of the choice of the lower scale \i. Secondly, at the level 
of the next-to-leading logarithm, we have 



B t {x,rj) =B t (x,l) + B QCD (v), 
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With 1] = 1.75 and rrit = 180 GeV obtained by averaging the data of CDF and DO groups 



15| , we have the values of the coefficient functions as follows 
A t = 2.09, B t 



-0.126, 
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b W r 2 



0.153. 
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Recent measurements of B — > by the CLEO group result in a constraint on the 
magnitude of F|, which is close to the value of the standard model. Likewise, the dilepton 
decay B — > X s l + l~ provides a unique test of the coefficients A and B which are absent in 
the radiative process. 

In addition to the top quark, the decay B — > X s l + l~ involves the charm quark and as 
well known, its total rate is actually dominated by the cc resonances. The contribution of 
these resonances such as J/ip and if)' to the effective Hamiltonian Eq. (1) may be taken into 
account by the usual form of vector-meson dominance [TE|-|rS|] 
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where My is the mass of the vector intermediate state and Ty its full width. The decay 
constant is defined such that fye^ = (Olcy^cl V(e)) and will be determined by the measured 
partial width for decays to lepton pairs fll9~] 
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with Q c = 2/3. In this work we treat the coupling for the neutral bscc four-quark operator 
0-2 = C\ + C2/3 as a phenomenological parameter and use the CLEO data | a-i \— 0.26 ± 0.03 
which is close to the value determined from a fit to the rate for the semi-inclusive 



process B — > X s J/ip [jnj]. The choice of a negative a 2 is in consistence with Breit-Wigner 
phase ip = [|22|| . However, consequences of assuming a different phase for a 2 can be derived 
without any difficulty in principle. Hence we can read off the vector resonance component 
for Eq. (2), 
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Moreover, there is an effect arising from quarks active at the \x scale, which contribute 
via penguin diagrams with insertion of four-quark operators. It is represented by 
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where 03 = C3/3 + C4, 04 = C3 + C4/3, a6 = C5 + C6/3, and the function comes from the 
one-loop matrix elements of the four-quark operators and is given by 
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Numerically, we have 03 = 0.0211, 04 = —0.0024, and a§ = 0.0028 for rj = 1.75. Comparing 
them with a 2 shows that the charm quark practically dominates the effect of these active 
quarks. Since the coupling of them to is of vector-type, the combination A — B remains 



unchanged after they are included. We remind the reader that the short-distance coefficients 
listed above are expressions in the NDR scheme. The scheme independence of the sum of 
these coefficients has been explicitly demonstrated by the authors of Ref . [[12] . 



III. FORWARD-BACKWARD ASYMMETRY IN DILEPTON B DECAYS 



A. B XJ+l- 



The differential forward-backward asymmetry of the l + production that we are to com- 
pute is defined by 

dT FB (q 2 ) = £ dT(cos6i) - J° dT(cos9i), 

in which 7r — Q\ is the polar angle of the l + with respect to the direction of motion of the 
decaying meson in the frame. For the inclusive process B — > X s l + l~ induced by the 
Hamiltonian in Eq. (1), it takes the form 
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which 



quarks respectively. The contribution of pure photonic penguin diagrams, 
governs the rate of radiative processes B — > X s 7, cancels out. Thus the measurement of the 
asymmetry provides entirely independent tests of the standard model. In Eq. (9) we have set 
the limit of mi — because there is no singularity at q 2 = 0; so there is no mixture between 
left- and right-handed leptons and the difference between electron and muon channels is 
negligible. As noticed by the authors of Ref. f23| , the forward-backward charge asymmetry 
is proportional to the coefficient combination (A — B) and may be large for the very heavy 
top quark. Recalling that (A — B) receives neither the QCD renormalization nor the cc 
effect, we anticipate that the asymmetry will become an effective measure for testing the 
high energy scale physics. To study its properties in detail, let us present a factorized form 
for the asymmetry {integrated over the interval of [0, q 2 ]) 
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Here I t stands for the contributions of loops containing the virtual top quark along with 
QCD corrections, for which we find 
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with s = q 2 jm\ and r s = m 2 /m 2 . In Fig. 1 we plot the integrated asymmetries arising from 
the photonic penguin and the Z penguin plus W box diagrams. They are rescaled in terms 
of the inclusive semileptonic decay width of the B meson 
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with fb c numerically equal to 0.39, in which we have included the one-loop QCD correction 
to the semileptonic B decay (all plots in this paper are rescaled in this way). For low 
s, photonic penguin diagrams dominate It(q 2 ), and an opposite sign of F\ to A t — B t in 
the standard model implies a positive asymmetry. As s increases, contributions of the Z 
penguin and W box diagrams emerge and cancel with that of photonic penguin diagrams. 
Consequently, the integrated asymmetry arising purely from virtual top quarks turns out to 
be zero at 

s = 1 
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For m t = 180 ± 12 GeV and r] = 1.75 we have s = 0.365 =f 0.010 



and s t — 

A t + B t 

which, by chance, is slightly below the J/ip peak point 0.400). The range of this zero 
point due to the error bar of the m t is expected to be small, considering that (A t + B t ) 
and F\ both depend rather weakly on the top mass. We find that this zero point varies 



from 0.356 =F 0.008 to 0.381 =F 0.012 for the above top quark masses when rj ranges from 
1.50 to 2.00, which incorporates uncertainties due to QCD scale parameter as well as the 
renormalization subtraction point. For the top component of the short-distance coefficient, 
the difference between NDR and HV schemes amounts to just about 5% ]13j. Thus we will 
see a slightly different value of the So in the HV scheme. Toward the high s region the Z 
penguin along with W box diagram takes over in magnitude and we are left with a surplus 
of the number of l~ scattered in the forward hemisphere in the / + /~ rest frame. 

Other components of the pure short-distance effect, coming from active quarks, are 
represented by 
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where the explicit form for IAz) is given in the Appendix and z = \ 1 -. Figure 2 

V s 

summarizes asymmetries coming from the short-distance effects alone. The nonresonant cc 
component dominating active quark effects becomes substantial above the threshold 2m c . 
Unfortunately, the total short- distance contribution suffers from a cancellation between the 
top and cc components, which reduces strongly the surplus of l~ of high q 2 in the forward 
hemisphere. With the same reasoning, we expect a vanishing net result of the short- distance 
effect as a whole somewhere in the middle region of the phase space. This sort of zero points 
(denoted by so) is identical in the NDR and HV schemes. To be precise, we evaluate the 
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position of this point for m t = 180 ± 12 GeV and r] = 1.75. It turns out to be s = 
0.511 =F 0.025, which comes closer to the ip' resonance than the J/ip. We find a weak rj 
dependence of this zero point and obtain a derivation less than 4% for rj varying from 1.50 
to 2.00 once again. 

Now let us consider resonant effects arising from vector intermediate states. We have 



with ry = My /mf and oy = T v /M v . The resonant component iy(<? 2 ) reaches its maximum 
value at s = ry. Because of the dependence on the real part of the resonant propagator, as 
s passes through r v , the differential asymmetry changes sign. There is , therefore, a partial 
cancellation of the Iv{.<f) and instead of a plateau we will find a decreasing integrated 
asymmetry for s across the mass of the vector resonance. 

The integrated asymmetry of the / + production in the inclusive process B — > X s l + l~ is 
plotted in Fig. 3 and includes the components of J/ip and ip' resonances. It is characterized 
by prominent peaks at q 2 = M 2 ^ and M 2 , respectively in the resonant region, as expected. 
The contribution of the ip' resonance is around one third that of the J/ip. In the J/ip resonant 
region, the top quark component of the asymmetry is actually small, as explained previously. 
It means that the product a 2 (A — B) dominates, as far as the coefficients are concerned. 
Unlike the decay rate, the (A — B) is involved here linearly. Alternatively stated, the short- 
distance effects as a whole disappears at s , which is near the ip'. The asymmetry, therefore, 
appears in the shape of the real part of the resonant propagator. Experimentally, a 2 has 
been reasonably well extracted from the CLEO data and hopefully more precise data will 
come from the planned B factories. As a result, the prediction for asymmetries near these 
zero points (or resonances) can be made in a way insensitive to the QCD renormalization 
parameter 77. Numerically, as r] ranges from 1.50 to 2.00, TpB(Mj^), for example, varies 
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merely by 6.1%. In our opinion, measurements of these asymmetries will be an effective way 
to unambiguously test the combination A — B arising purely from short- distance physics. 
Meanwhile, as A — B depends strongly on the top mass in the standard model, it may be 
an independent way to confirm the value of rrit determined in different ways. 



B. B -» K*l + r 

We now turn our attention to the exclusive channel B — > K* . The effective quark 
current sY^ B 'b in question has two different structures; the parametrization for the matrix 
element of V — A currents in terms of invariant form factors is 



(p>|*y"(l-7 5 )&|P>= [a + (q 2 )(P + p), + a„(q 2 )(P-p),]P 
In analogy to this we have for the magnetic-moment operator 
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+f(q 2 )<i>; + ig(q 2 )e^ aP r u P a p , 

along with a condition of current conservation, 

(M 2 - m 2 )a+ + g 2 a„ + / = 0. 

Here M (P M ) and m (p^) are masses (momenta) of B and K* mesons, respectively, </> M the 
polarization vector of the K* meson (satisfying (j)^ = 0) and q = P — p the momentum 
transfer into the dilepton. 

The differential forward-backward asymmetry of B — > K* reads || 
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~dq^ = U^\^2) [atts 2 w J TP 

where necessarily now + m 2 — q 2 ) 2 — AM 2 m 2 . The helicity amplitudes ap- 

pearing here are defined as 

Hi {K) = A(B)h ± + m b s 2 w Ffh ±} (18) 
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where h± = f ± XMg |24| and h± = f ± XMg. Assuming the B meson contains an on-shell 
b quark of velocity v, one finds 
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;(M-VPTX) h±. 



(19) 



(see Ref. 0,|J for detailed discussion.) With these relations we may rewrite the differential 
asymmetry of Eq. (16) as 
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Certainly, further analysis requires the knowledge of hadronic form factors, but it is 
hard to calculate them directly from the first principles of QCD. In Ref. |5],|6j we have used 
heavy flavor symmetry to relate the form factors for B — > K*l + l~ decays with those for 



D — > K* l + v. In the leading order of heavy quark effective theory, form factors scale as |25 
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la this equation form factors are evaluated at the same v ■ p, namely that / and g at 
q 2 = M 2 + m 2 — 2Mv ■ p are related to fo and go at q 2 D = M|, + m 2 — 2M D v ■ p. Thus 
the whole physical range of q 2 D = [0, (Md — m) 2 ] in D — > covers g 2 = [q 2 , (M — m) 2 } 
in B — >• if*, where = (4.07 GeV) 2 . With the data available for — > K* we are led to 
a reliable estimate of the asymmetry in the region above the ip' resonance, on the basis of 
momentum-to-momentum correspondence of form factors 0. However, extrapolating the 
data of D decays to lower q definitely needs the knowledge of the g 2 -dependence for form 
factors. In this work we will employ the pole dominance approximation, 
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h(q 2 ); (h = f,g). 
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Here M' is the pole mass of the current involving b and s quarks and we make no distinction 
between the vector and axial vector masses. While this form seems reasonable for the small 
recoil region of B — > K* decays, as the end-point is close to the pole, we hope it helps to 
figure out the properties of forward-backward asymmetries in the resonant region or below. 
Regarding poles of degree n, we make use of two schemes; one assumes the monopole form 
for both / and g like the BSW model |[26|| , and another follows the suggestion of Korner and 
Schuler ||27|| , namely n = 1 for /, but a dipole form of n = 2 for g (The latter is referred to 
as the dipole scheme in this paper despite / actually being of monopole form.). Meanwhile, 
Eq. (20) will be used to obtain form factors at the kinematical point go; so with form 
factors Ai(0) and V(0) (they are linear combinations of fo and gr>) of D — > K*(892) l + v 
measured by E691, E687 and CLEO groups 0], one ends up with /r>(0) = —1.43 GeV and 
g D (0) = 0.695 GeV" 1 , which translate to f{ql) = -2.64 GeV and g{ql) = 0.455 GeV" 1 . 

Before doing practical calculations, we perform a simple analysis of integral asymmetries, 
by using the monopole form factor along with the limit of M = M' and m 2 /M 2 < 1. It 
ought to be appropriate for the region not close to the pole, say around the J/ip mass or 
below. In this case, the counterpart of Eq. (11) in the K* mode, with necessarily s = q 2 /M 2 
for this moment, has the form 
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We have defined that fo = /(0) and go = g(0). Once again, competition between photon 



penguin and Z penguin plus W box diagrams causes a zero point of I t 
,(2 ~st) fo > 
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Clearly, as s t < 1 in the standard model, the condition that fo is of opposite sign to go 
guarantees a real value of this root. This condition is satisfied physically because the helicity 
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relation | h+ \<\ h~ | has been well confirmed in semileptonic decays p4j| , namely that the 



V — A structure of the quark current is manifested at the hadron level. Numerically, with 
/(gg) and g{q^) given earlier the zero point turns out to be 0.345=f0.07 for m t = 180±12 GeV 
and 1] = 1.75. It coincides with the J/ip peak M 2 ^/M 2 = 0.344. 
Meanwhile, the correspondence to Eq. (12) reads 



M 5 \ U t \ 2 ( G, x 2 



TV) = =^ (T2 f ° g ° 



(oa + a 4 + fle) & K *\m 2 c /M 2 , s) - ^ $^(0, s) 

2 (25) 

03 + 04 + CLq fo(K*)(-\ \ a 4 + a 6 2 



with 



$ (Ar,) (r 9 , s) = Ar q s - i(lnr, - ^)s 2 - 16r 9 2 [2/ 3 (^) - / 2 (z) - (26) 

In analogy to the inclusive process, adding the active quark effect to the top one ends 
up with a cancellation of the asymmetry once again. We find a complete cancellation at 
So = 0.451 =F 0.012 for the same top mass and rj as before. In addition we work out the 
vector resonant contribution which is 

4r V) - g§) 2 + H , (^) 

where now ry = M 2 /M 2 . We should remark that the above expressions illustrate the 
features of the asymmetries in 5 — > K*l + l~ , viz., the peaks in the resonant contribution, 
the cancellation between Z penguin and plus W box diagrams, and the plateau in the 
nonresonant cc effect above the threshold 2m c . In particular, we found the position of the 
zero point for the short- distance effect is in good agreement with what we obtained via the 
full expression. 

Nevertheless, if we still ignored the K* mass in the region near the endpoint, the phase 
space would approach q 2 = M 2 rather than (M — m) 2 . While it may make a slight difference 
in the inclusive process (where m is replaced by m s ), this small change in the upper boundary 
of phase space could lead to a sizeable overestimate of observable quantities such as the decay 
rate and asymmetry in exclusive channels. This is because the pole in the form factors 
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dominates the distribution at the endpoint. The shift of the pole mass to M would further 
increase the estimate. Thus when the whole range of q 2 in the decay B — > K* is considered, 
we apply the full expression of Eq. (19) (using the experimental values for the masses m 
and M'). In this way, the integrated asymmetry for B — > K*l + l~ is plotted in Fig. 4 for 
monopole and dipole form factors. We have checked that there is a variation of the zero 
point of short-distance effects due to the choice of form factors. With the same m t and r\ 
as before, there are s = 0.338 =F 0.007 (the top quark effect alone) and s = 0.511 =F 0.017 
(the short-distance effect as a whole) for the monopole scheme, as well as 0.291 =F 0.006 
and 0.444 =p 0.013 for the dipole one. The contribution of J/ip along with the cc effect still 
dominates its resonant region. By comparison, we find the two schemes of form factors 
give similar results of the integrated asymmetry above the ip' . This is expected when both 
schemes use the data of semileptonic D decays as input. But asymmetries, integrated over 
the interval between the J / ip and ip' and the region below the J / ip with the monopole form, 
are 1.7 times as large as that with the dipole one. This factor becomes even larger when q 2 is 
near to the origin. For the monopole form, a fraction of 25% or so of the asymmetry occurs 
in the ground state and this is far from saturating the inclusive process. So it is anticipated 
that higher if-resonances may play an important role in measuring asymmetry as well. 



IV. DISCUSSION AND CONCLUSION 

As we have seen in Fig. 1, the magnetic- moment operator becomes important for photons 
approaching the mass-shell. The positive asymmetry arising from it is the consequence of 
the coefficient F\ being negative relative to A — B in the standard model. Thus measurement 
in this region will tell us the sign of Fj, as well as confirming its bound given by the decay 
b — > 57. 

On the other hand, the relative phase between F\ and A t + B t impacts dramatically on 
the b — > sl + l~ . For example, a potential magnetic-moment operator of a similar size but 
opposite sign to that in the standard model, such as those allowed in supersymmetry models 
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29fl , may increase the forward-backward asymmetry of the whole kinematical region by a 
factor of 4.0 ~ 5.0 and a factor 2.0 for that in the low q region with a cut-off at Mj/^, but 
reversing the sign of the asymmetry. In fact, if this relative phase is negative it shows us 
evidence of new physics, whereas if it is positive it leads to the sensitivity of the asymmetry 
with respect to the J /if} and the cc continuum in the middle region. We emphasize that 
the relative phase between F\ and A t + B t manifests itself as well in the invariant mass 
spectrum of dileptons. When left-handed leptons do not mix up with right-handed ones in 
the massless lepton limit, the only interference occurs between quarks outgoing from the four- 
fermion vertex and that from the penguin, yielding a term proportional to Re[(A t + B t )F^\ 
in the dilepton spectrum. As it is likely that the first experiment will be carried out for 
this spectrum, we will be able to see whether there is cancellation between the top quark 
components of the asymmetry in a model independent way. 

Toward high q 2 , coefficients A t and B t become important as they are enhanced by an 
extra q 2 . The contributions of Z boson penguin along with W box diagrams in this region 
surpass that of the photonic penguin diagram for a large top quark mass such as 174 GeV 
(see Eq. (6)). While A — B (and Ft/) may be determined by the method described before, 
this asymmetry will help to extract the value of the combination A t + B t experimentally. 
When B t always occurs in combinations with the much larger A t , any appreciable difference 
between A t + B t and A t — B t is beyond the standard model. 

As for the exclusive decay B — > K*l + l~, the difference between the two assumptions 

about form factors becomes evident for low q. When applied to the rare radiative decay 

B — > K*7, the monopole form factors give a theoretical value of 17% for the ratio of the 

Y(B -> K*7) 

exclusive to inclusive rate, defined as R = — -. This is favored by the preliminary 

T(B -> X s >y) 

CLEO data |IJ, while the dipole scheme result of 5.7% is not. Although these two sets of 
form factors are similar to each other near the endpoint, we tend to favor the prediction for 
charge asymmetries obtained by using the monopole scheme, as far as the whole kinematical 
region is concerned. 

Finally we want to point out that possible experimental cuts in the invariant mass of 
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lepton-antilepton pairs do not necessarily mean loss of statistics as far as asymmetry is con- 
cerned. For instance, a cut-off at the J/ip mass from which we integrate to the endpoint 
<?max leads to an asymmetry larger in magnitude than that over the whole kinematical re- 
gion. In fact, this interval gives the maxinum integrated asymmtry because of the resonant 
enhancement. Our prediction for it, using m t = 180 GeV, is — 1.54 x 10 _5 r s l . in the inclusive 
process, which is accessible in the upcoming B factories. 

In conclusion we have presented a comprehensive analysis of forward-backward asymme- 
try of l + production in the decay B — > X s l + l~ and B — > K*l + l~. We have attempted to 
specify certain regions in phase space of the B decay in which the asymmetry is sensitive to 
individual short-distance coefficients. In particular, we note the potential importance of the 
J /if), since the integrated asymmetry near its resonant region is largely independent of the 
QCD renormalization. Given the strong dependence of A — B on the top quark mass, mea- 
surements of this asymmetry will not only provide an effective test of the standard model 
but also lead to information about m t . As well, we believe that the shape of the plots in Figs. 
1 to 4 represent well the expectation for the forward-backward asymmetry in the framework 
of the standard model. Therefore if future experiments observe any dramatic deviation from 
them, they can be regarded as the evidence of new physics beyond the standard model. 
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APPENDIX: INTEGRAL FORMULAE 



When the effect of the one-loop matrix element is studied in the main text, we have used 
the notation Ik(z) for the integral defined as 



I k (z) = Re 



du 



,1 — U 



2\k 



. 1 +U . 

in in 

1 — u 



(fc = l,2,...). 



It holds in both ranges above and below the threshold 2m q . In the latter case where u is 
imaginary, we apply the analytic continuation 



arctan 



u 



, 1 + u . 

in nx 

l — u 



The lower limit z corresponds to q 2 = 0. It is not difficult to work out Ik for k — 1 



h{z) = ^Re 



. 1 + u ( 1 , 1 + u . 

In - In in 

1 — u \2 1 - u 



zo 



We would like to present in this appendix a general procedure to carry out the integral for 
arbitrary k. To this end, it is convenient to consider an auxiliary function 



T(z; a) = Re 



du 



z (a — v?) 



, 1 + u . 
m m 



u 



(a>0), 



a=l • 



and thus the integral in question may be rewritten as 

(-l)(*-^(*-Dj(z;q) 
(jfe-1)! Sa^ 1 ) 

Hence it is straightforward to derive the Ik(z) for any k larger than one with the explicit 
form of I(z; a) in terms of a series 

1 1 



+ 



1 



(1 
1 



[l + uy 



u) 1 



, 1 +M . s 

(In m) 



u 
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in which Jj^ is given by the recurrence 

du 1 1 



la 



(1±u)'(1tu) 2 
where we know the starting values 



=F 



(Z-l)(l±«) 



1 1 



it 
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Figure captions: 

Figure 1 The top quark components of forward-backward charge asymmetry in the 
decay B — > X s / + Z~, integrated over the region of [0, q 2 }. The dashed line corresponds to 
the photon penguin diagram, while the solid to Z penguin plus W box diagrams. We take 
the values m t = 174 GeV, i] = 1.75, and m b = 4.9 GeV and the limit of m 2 s jm\ — > 0. Both 
plots are rescaled by the semileptonic decay width of the B meson, multiplied by the factor 
of 10 5 . 

Figure 2 Integrated asymmetries due to short-distance effects in B — > X s . The 
dashed line corresponds to the active qq contribution dominated by the cc, and the solid 
to the top quark one, that is the sum of plots in Fig. 1. We use m c = 1.5 GeV and the 
remaining parameters and the scale are the same as in Fig. 1. 

Figure 3 The same as in Figs. 1 and 2 for the sum (dashed line) of all effects (see Eq. 
(10) in the main text), including resonant contributions (solid line). 

Figure 4 Integrated forward-backward asymmetries in the decay B — > K* The 
dashed line is associated with the monopole scheme of form factors, while the solid line goes 
with the dipole one. We take M = 5.28 GeV, m = 0.892 GeV, and M' = 5.38 GeV. The 
remaining parameters and the scale are the same as in previous figures. 
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